A LOCAL-TIME CORRESPONDENCE FOR 
STOCHASTIC PARTIAL DIFFERENTIAL EQUATIONS 
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Abstract. It is frequently the case that a white-noise-driven parabolic and/or hyperbolic 
stochastic partial differential equation (SPDE) can have random-field solutions only in spa- 
tial dimension one. Here we show that in many cases, where the "spatial operator" is the 
L 2 -generator of a Levy process X, a linear SPDE has a random- field solution if and only if 
the symmetrization of X possesses local times. This result gives a probabilistic reason for 
the lack of existence of random-field solutions in dimensions strictly bigger than one. 

In addition, we prove that the solution to the SPDE is [Holder] continuous in its spatial 
variable if and only if the said local time is [Holder] continuous in its spatial variable. We also 
produce examples where the random-field solution exists, but is almost surely unbounded in 
every open subset of space-time. Our results are based on first establishing a quasi-isometry 
between the linear L 2 -space of the weak solutions of a family of linear SPDEs, on one hand, 
and the Dirichlet space generated by the symmetrization of A, on the other hand. 

We study mainly linear equations in order to present the local-time correspondence at a 
modest technical level. However, some of our work has consequences for nonlinear SPDEs 
as well. We demonstrate this assertion by studying a family of parabolic SPDEs that have 
additive nonlinearities. For those equations we prove that if the linearized problem has 
a random-field solution, then so does the nonlinear SPDE. Moreover, the solution to the 
linearized equation is [Holder] continuous if and only if the solution to the nonlinear equation 
is. And the solutions are bounded and unbounded together as well. Finally, we prove that 
in the cases that the solutions are unbounded, they almost surely blow up at exactly the 
same points. 
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1. Introduction 

We consider the stochastic heat equation inspired by the fundamental works of Pardoux 
(I1975al 11975b! fT972l . Krylov and Rozovskii (I1979al I1979H fT977j) . and Funaki ffl983l . Let 
w denote space-time white noise, t is nonnegative, x is in R^, and the Laplacian acts on the 
x variable. Then, we have: 



(1.1) 



d t H(t,x) = (AH)(t,x) + w(t,x), 
H(0,x) = 0. 

Let us consider also the stochastic wave equation of Cabana (119701) : 

' d tt W(t,x) = (AW)(t,x)+w(t,x), 

W(0,x) = dtW(0,x) = 0. 

One of the common features of ( 11. lj) and ( 11.21) is that they suffer from a curse of dimen- 
sionality. Namely, these equations can have random-field solutions only in dimension one. 
Moreover, this curse of dimensionality appears to extend beyond the linear parabolic setting 



;i-2) 
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of (11.11) . or the linear hyperbolic setting of (jl.2p . For instance, see Perkins (20021 Corollary 
III. 4. 3) for an example from superprocesses, and Walsh (119861 Chapter 9) for an example 
from statistical mechanics. 

One can informally ascribe this curse of dimensionality to the "fact" that while the Lapla- 
cian smooths, white noise roughens. In one dimension, the roughening effect of white noise 
turns out to be small relative to the smoothing properties of the Laplacian, and we thus 
have a random-field solution. However, in dimensions greater than one white noise is much 
too rough, and the Laplacian cannot smooth the solution enough to yield a random field. 

Dalang and Frangos (119980 were able to construct a first fully-rigorous explanation of the 
curse of dimensionality. They do so by first replacing white noise by a Gaussian noise that is 
white in time and colored in space. And then they describe precisely the roughening effect 
of the noise on the solution, viewed as a random generalized function. See also Brzezniak 
and van Neerven (12003]) . Dalang and Mueller ([2003]) . Millet and Sanz-Sole (1T9991) . Peszat 
(12002]) . and Peszat and Zabczyk (120001) . More recently, Dalang and Sanz-Sole (12005]) study 
fully nonlinear stochastic wave equations driven by noises that are white in time and colored 
in space, and operators that are arbitrary powers of the Laplacian. 

In this article we present a different explanation of this phenomenon. Our approach 
is to describe accurately the smoothing effect of the Laplacian in the presence of white 
noise. Whereas the answer of Dalang and Frangos (1998) is analytic, ours is probabilistic. 
For instance, we will see soon that ( 11. ip and ( 11.21) have solutions only in dimension one 
because d- dimensional Brownian motion has local times only in dimension one [Theorem 
12.1] . Similarly, when d = 1, the solution to (11.11) [and/or (11.21) ] is continuous in x because 
the local time of one-dimensional Brownian motion is continuous in its spatial variable. 

The methods that we employ also give us a local-time paradigm that makes precise the 
claim that the stochastic PDEs (11.11) and (11.21) "have random-field solutions in dimension 
d = 2 — e for all e G (0 , 2)." See Example 17.41 below, where we introduce a family of SPDEs 
on fractals. 

An outline of this paper follows: In §|2]we describe a suitable generalization of the stochastic 
PDEs (11.11) and ( 11.21) that have sharp local-time correspondences. This section contains the 
main results of the paper. The existence theorem of $2] is proved in ^J3j Section H] contains 
the proof of the necessary and sufficient local-time condition for continuity of our SPDEs in 
their space variables. Although we are not aware of any interesting connections between local 
times and temporal regularity of the solutions of SPDEs, we have included $5] that contains 
a sharp analytic condition for temporal continuity of the SPDEs of §0 We also produce 
examples of SPDEs that have random-field solutions which are almost surely unbounded in 
every open space-time set [Example 15.5] . 
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Section [6] discusses issues of Holder continuity in either space, or time, variable. In £J7] we 
establish a very general (but somewhat weak) connection between Markov processes, their 
local times, and solutions to various linear SPDEs. The material of that section is strongly 
motivated by the recent article of Da Prato (120070 who studies Kolmogorov SPDEs that 
are not unlike those studied here, but also have multiplicative nonlinearities. We go on to 
produce examples where one can make sense of a random-field solution to (11.11) in dimension 
2 — e for all e G (0 , 2). These solutions in fact turn out to be jointly Holder continuous, but 
we will not dwell on that here. Finally, in £j8]we discuss a parabolic version of the SPDEs of 
£j2]that have additive nonlinearities; we prove the existence of solutions and describe exactly 
when and where these solutions blow up. 

The chief aim of this paper is to point out various interesting and deep connections be- 
tween the local-time theory of Markov processes and families of stochastic partial differential 
equations. In all cases, we have strived to study the simplest SPDEs that best highlight these 
connections. But it would also be interesting to study much more general equations. 

There are other connections between local times of Markov processes and Gaussian pro- 
cesses that appear to be different from those presented here. For a sampler of those isomor- 
phism theorems see Brydges, Frohlich, and Spencer (11982[) . Dynkin ( 119841) . and Eisenbaum 
(I1995p . Marcus and Rosen (120060 contains an excellent and complete account of this theory. 
Diaconis and Evans (120021) have introduced yet a different isomorphism theorem. 

Finally, we conclude by mentioning what we mean by "local times," as there are many 
[slightly] different versions in the literature. Given a stochastic process Y := {Y t } t > on R d , 
consider the occupation measure [s], 



We can identify each Z(t , •) with a measure in the usual way. Then, we say that Y has local 
times when Z(t , dx) <C dx for all t. The local times of Y are themselves defined by Zit , x) : = 
Z(t , dx)/dx. It follows that if Y has local times, then Z(t , (p) = j Rd Zit , x)ip(x) dx a.s. for 
every t > and all measurable functions (p : R d — ► R + . And the converse holds also. 



In order to describe when ( 11.11) and ( 11. 2ft have random-field solutions, and why, we study 
more general equations. 

Let Jzf denote the generator of a <i-dimensional Levy process X := {X t }t>o with charac- 
teristic exponent We can normalize things so that Eexp(i£ ■ X t ) = exp(— and 



(1.3) 




tp(Y s ) ds for all t > and measurable cp : H d — ► R + . 



2. More general equations 
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consider Jzf as an L 2 -generator with domain 



(2-3) 



(2.1) DomJ*f := j/ G X/(R«) : jf |/(OH*(OI < °° 
As usual, / denotes the Fourier transform of /; we opt for the normalization 

(2.2) /(£) := I e^ x f(x) dx for all / G L 1 (R d ). 

[The general L 2 -theory of Markov processes is described in great depth in Fukushima, 
Oshima, and Takeda (119940 in the symmetric case, and Ma and Rockner (JT992) for the 
general case.] 

In this way, we can — and will — view Jzf as a generalized convolution operator with Fourier 
multiplier Jzf(£) := -*(£)■ 

We consider two families of stochastic partial differential equations. The first is the sto- 
chastic heat equation for Jzf : 

d t H(t,x) = (3fH)(t,x)+w(jt,x), 
H{0,x) = 0, 

where x ranges over R d and t over R + := [0 , oo), and the operator Jzf acts on the x variable. 
We also consider hyperbolic SPDEs of the following wave type: 

d tt W(t,x) = (S?W)(t,x)+w(t,x), 

W(0,x) = d t W{0,x) = 0. 

Recall that X := {X t } t > is a Levy process whose generator is Jzf. Let X' denote an 
independent copy of X and define its symmetrization, a la Paul Levy, by 

(2.5) X t := X t - X' t for all t > 0. 

It is a standard fact that X is a symmetric Levy process with characteristic exponent 2Re\l/. 

The following is our first main result, though the precise meaning of its terminology is yet 
to be explained. 

Theorem 2.1. The stochastic heat equation (12. 3p has random-field solutions if and only if 
the symmetric Levy process X has local times. The same is true for the stochastic wave 
equation (j2.4j) . provided that the process X is itself symmetric. 

Remark 2.2. Brzezniak and van Neerven (I2U03I) consider parabolic equations of the type 
( 12.31) . where Jzf is a pseudo-differential operator with a symbol that is bounded below in 
our language [bounded above in theirs]. If we apply their theory with our const ant- symbol 
operator Jzf, then their bounded-below condition is equivalent to the symmetry of the Levy 
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process X, and their main result is equivalent to the parabolic portion of Theorem 12. 11 under 
the added condition that X is symmetric. □ 



It is well-known that when d>2, Levy processes in R d do not have local times (Hawkes, 1986) 



It follows from this that neither (12.31) nor (12.4ft [under a symmetry assumption on X] can 
ever have random-field solutions in dimension greater than one. But it is possible that there 
are no random-field solutions even in dimension one. Here is one such example; many others 
exist. 

Example 2.3. Suppose X is a strictly stable process in R with stability index a G (0,2]. 
It is very well known that the symmetric Levy process X has local times if and only if 
a > 1; see, for example, Hawkes (11986). If X is itself symmetric, then Jzf = — (— A) a / 2 is the 



a-dimensional fractional Laplacian (Stein, 1970 Chapter V, §1.1), and also the stochastic 



wave equation (12 Ah has a random-field solution if and only if a > 1. □ 

The local-time correspondence of Theorem 12.11 is not a mere accident. In fact, the next 
two theorems suggest far deeper connections between the solutions to the linear SPDEs of 
this paper and the theory of local times of Markov processes. We emphasize that the next 
two theorems assume the existence of a random-field solution to one of the stochastic PDEs 
(12. 3p and/or (12. 4p . Therefore, they are inherently one dimensional statements. 

Theorem 2.4. Assume that d = 1 and the stochastic heat equation (12.31) has a random-field 
solution {H(t , x)}t>o,x£R- Then, the following are equivalent: 

(1) There exists t > such that x i— > H(t ,x) is a.s. continuous. 

(2) For all t > 0, x h- > H(t ,x) is a.s. continuous. 

(3) The local times of X are a.s. continuous in their spatial variable. 

The same equivalence is true for the solution W to the stochastic wave equation (12 .4p . pro- 
vided that the process X is itself symmetric. 

Theorem 2.5. Assume that d = 1 and the stochastic heat equation (12.31) has a random-field 
solution {H (t , x)} t >o,xe~R- Then, the following are equivalent: 

(1) There exists t > such that x i— > H(t ,x) is a.s. Holder continuous. 

(2) For all t > 0, x i— > H(t , x) is a.s. Holder continuous. 

(3) The local times of X are a.s. Holder continuous in their spatial variable. 

If the process X is itself symmetric, then the preceding conditions are also equivalent to 
the Holder continuity of the solution W to the stochastic wave equation ( 12.41) in the spatial 
variable. Finally, the critical Holder indices of x ^ H(t,x), x i— > W(t,x) and that of the 
local times of X are the same. 
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Blumenthal and Getoor (11961) have introduced several "indices" that describe various 
properties of a Levy process. We recall below their lower index f3": 

2.6 p := limmt — — ; — = sup < a > : hm — — = oo 

Theorem 2.6. If /3" > d, then the stochastic heat equation (12. 3p has a random-field solution 
that is jointly Holder continuous. The critical Holder index is < {(3" — d)/2 for the space 
variable and < (/3" — d)/2/3" for the time variable. Furthermore, if X is symmetric, then the 
same assertions hold for the solution to the stochastic wave equation i\2A\\ . 

Remark 2.7. If X is in the domain of attraction of Brownian motion on R, then f3" = 2, 
and the critical temporal and spatial index bounds of Theorem 12.61 are respectively 1/2 and 
1/4. These numbers are well known to be the optimal Holder indices. In fact, the relatively 
simple Holder-index bounds of Theorem 12.61 are frequently sharp; see Example 15.41 □ 



It is well known that < f3" < 2 ( Blumenthal and Getoor, 1961 Theorem 5.1). Thus 



the preceding is inherently a one-dimensional result; this is in agreement with the result of 
Theorem 12. 11 Alternatively, one can combine Theorems 12.11 and 12.61 to construct a proof of 
the fact that (3" < 2. [But of course the original proof of Blumenthal and Getoor is simpler.] 

Several things need to be made clear here; the first being the meaning of a "solution." 
With this aim in mind, we treat the two equations separately. 

We make precise sense of the stochastic heat equation (12.31) in much the same man- 
ner as Walsh ffTM : see also Brzezniak and van Neerven (120031) . Dalang fTl999l I20011 . 
and Da Prato (2007). It is well known that it is much harder to give precise meaning 
to stochastic hyperbolic equations of the wave type (12 .4p . even when Jz? is the Laplacian 



(Dalang and Leveque, 2004a , 2004b Dalang and Mueller, 2003 Quer-Sardanyons and Sanz-Sole, 2003 



Mueller, 1997| |Dalang, 1999[ [Mueller, 1993| ). Thus, as part of the present work, we introduce 
a simple and direct method that makes rigorous sense of (12.41) and other linear SPDEs of 
this type. We believe this method to be of some independent interest. 

2.1. The parabolic case. In order to describe the meaning of (12.31) we need to first intro- 
duce some notation. 

Let {Pt}t>o denote the semigroup of the driving Levy process X; that is, (P t f)(x) : = 
Ef(x + X t ) for all bounded Borel-measurable functions / : R d — > R [say], all x G R d , and 
all t > 0. [As usual, X := 0.] Formally speaking, P t = exp(tJzf). 

The dual semigroup is denoted by P*, so that (P 4 */)(x) = E/(x — X t ). It is easy to see 
that P and P* are adjoint in L 2 (R a! ) in the sense that 

(2.7) (P t f,g) = (f,P t *g) for all /, g e L 2 (K d ). 



8 FOONDUN, KHOSHNEVISAN, AND NUALART 

Needless to say, (• , ■) denotes the usual Hilbertian inner product on L 2 (R d ). 

Let y(R d ) denote the class of all rapidly decreasing test functions on R d , and recall that: 
(i) y(R d ) C Dom(^) n C°°(R d ) C L 2 (R d ) n C°°(R d ); and (ii) for all if, ip G y(R d ) } 

(2.8) Km( Pt+hLp - PtLp ,A=(J?P tV ^). 



HO \ h 

Thus, v(t , x) := (P t <f)(x) solves the Kolmogorov equation d t v(t , x) = {S£v){t , x), subject to 
the initial condition that v(0 ,x) = if(x). This identifies the Green's function for dt — Jzf = 
0. Hence, we can adapt the Green-function method of Walsh (1986, Chapter 3), without 
any great difficulties, to deduce that a weak solution to (12. 31) is the Gaussian random field 
{H{t , if); t > 0, if G y(R d )}, where 



(2.9) H(t,<p):= f f (P;_ sV )(y)w(dyds). 

JO JR d 



This is defined simply ClS db Wiener integral. 

Proposition 2.8. The Gaussian random field {H(t,ip); t > 0, if G J^(R d )} is well defined. 
Moreover, the process if i— > iY(t , <£>) zs a.s. linear for each t > 0. 

Proof. On one hand, the Wiener isometry tells us that 

(2.10) E(|tf(t,^)| 2 ) = t \\P;_ s f\\ 2 LH K*)ds. 

Jo 

On the other hand, it is known that each P* is a contraction on L 2 (R d ). Indeed, it is not 
hard to check directly that if id denotes the Lebesgue measure on R d , then the dual Levy 



process —X is ^-symmetric (Fukushima et al., 1994, pp. 27-28). Therefore, the asserted 



contraction property of P* follows from equation (1.4.13) of Fukushima, Oshima, and Takeda 
(T19M p. 28). It follows then that E(\H(t , f)\ 2 ) < t\\ip\\ 2 L2(Rdv and this is finite for all 
t > and if G y(R d ). This proves that u is a well-defined Gaussian random field indexed 
by R + x '(R d ). The proof of the remaining property follows the argument of Dalang 
(|1999j Section 4) quite closely, and is omitted. □ 

2.2. The nonrandom hyperbolic case. It has been known for some time that hyperbolic 
SPDEs tend to be harder to study, and even to define precisely, than their parabolic coun- 
terparts. See, for instance, Dalang and Sanz-Sole (120070 for the most recent work on the 
stochastic wave equation in dimension 3. 

In order to define what the stochastic wave equation (12.41) means precisely, we can try 
to mimic the original Green-function method of Walsh (119861) . But we quickly run into 
the technical problem of not being able to identify a suitable Green function (or even a 
measure) for the corresponding integral equation. In order to overcome this obstacle, one 



(2.11) 
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could proceed as in Dalang and Sanz-Sole (I2007p . but generalize the role of their fractional 
Laplacian. Instead, we opt for a more direct route that is particularly well suited for studying 
the SPDEs of the present type. 

In order to understand (12.41) better, we first consider the deterministic integro-differential 
equation, 

d tt u(t,x) = (Sfu)(t,x)+f(t,x), 
u(0 , x) = d t u(0 , x) — 0, 

where / : R + x R d — > R is a "nice" function, and the variables t and x range respectively 
over R + and R d . We can study this equation only under the following symmetry condition: 

(2.12) The process X is symmetric. 

Equivalently, we assume that \I/ is real and nonnegative. 

Recall that denotes the Fourier transform in the x variable, and apply it informally 
to (12. lip to deduce that it is equivalent to the following: For all t > and £ G R d , 

d tt u(t,o = -y(tWt,0 + f(t,0, 
w(o,o = $«(o,0 = o- 

This is an inhomogeneous second-order ordinary differential equation [in t] which can be 
solved explicitly, via Duhamel's principle, to produce the following "formula" : 

(2. 14) u(t , = -/== J* sin ( y/¥ffi (t - 8 )) f\s , ds. 

We invert the preceding — informally still — to obtain 

! f *sm(y/*(£j(t-sj) A 

(2.15) u(t,x) = ——, / ^ == J -e-^ x f(s,i)dsdi. 

We can multiply this by a nice function 99, then integrate [dx] to arrive at 



(2.13) 



1 r rtsm(y/*(g)(t 



(2.16) u(t , <p) := -— d / / ^ == (p{i) f(s , ds dt 

We may think of as this as the "weak/distributional solution" to (12.111) . 



2.3. The random hyperbolic case. We follow standard terminology and identify the 
white noise w with the iso-Gaussian process {iy(^)}/ ie z,2( R+xR d) as follows: 

p /»oo 

(2.17) w(h) := / / h(s , x) w(ds dx). 

Jn d Jo 
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Next, we define the Fourier transform w of white noise: 
(2.18) »■(») := Jggj = ^-L. ^ jT°S(.,{)«(<i.de), 

all the time remembering that "/i" refers to the Fourier transform of ft, in its spatial variable. 
Suppose h(s ,x) = ipi(s)ip 2 (x) for t > and x G R d , where y>i G L 2 (R + ) and ip 2 G L 2 (R d ). 
Then, it follows from the Wiener isometry and Plancherel's theorem that 

( 2 - 19 ) IK^IIi^xRd) = ii^iii2 (R+XRd) . 

Because L 2 (R + ) ®L 2 (Tl d ) is dense in L 2 (R + x R d ), this proves that w is defined continuously 
on all of L 2 (R + x R d ). Moreover, w corresponds to a white noise which is correlated with 
w, as described by the following formula: 

(2.20) E [w{h 1 )- : ^}=—^f J^h{s~0h 2 (s,0^ds, 

valid for all h 1 ,h 2 e L 2 (K + x R d ). 

In light of (12.161) . we define the weak solution to the stochastic wave equation (12.41) as the 
Wiener integral 



1 



sin(V^)(t- S ; 



(2-21) Wit^^JJ, ^ 

It is possible to verify that this method also works for the stochastic heat equation (12.31) . and 
that it produces a equivalent formulation of the Walsh solution (12. 9p . However, in the present 
setting, this method saves us from having to describe the existence [and some regularity] of 
the Green function for the integral equation (12.111) . 

Proposition 2.9. If the symmetry condition fl 2 . 1 2 [) holds, then the stochastic wave equation 
AMD has a weak solution W for all tp G ^(R d ). Moreover, {W(t,<p); t > 0, <p G ^(R d )} 
is a well-defined Gaussian random field, and <p i— > W(t , cp) is a.s. linear for all t > 0. 

Proof. We apply the Wiener isometry to obtain 

(2.22) e(| H -( ( ,^)I 2 ) = If — ^ J -mw^. 

Because | sin6/6\ < 1, 

E(\W(t jV )\ 2 ) < JL [ [\t-s) 2 M0\ 2 dsdt 
{ A7l > Jn d Jo 



(2.23 



2 



~^\\ l P\\L 2 (R d )i 
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thanks to Plancherel's theorem. It follows immediately from this that {W(t ,<p); t > 0, (p G 
S^(R d )} is a well-defined Gaussian random field. The remainder of the proposition is stan- 
dard. □ 



3. Existence of functions- valued solutions: Proof of Theorem 12. ll 

Let u := {u(t,ip); t > 0, <p G y(R d )} denote the weak solution to either one of (12.31) or 
(12. 4p . Our present goal is to extend uniquely the Gaussian random field u to a Gaussian 
random field indexed by R + x M, where M is a maximal subset of 2t(R d ) — the space of all 
Schwartz distributions on R d . Such an M exists thanks solely to functional- analytic facts: 
Define, temporarily, 

(3.1) d t (<p) := yjE(\u(t,(p)\ 2 ) for all ip G y(R d ) and t > 0. 

Then, the linearity of u in (p shows that (ip i— > d t (<p — ip) defines a metric for each t > 0. 
Let M t denote the completion of y(R, d ) in 3t(R d ) with respect to the metric induced by df, 
and define M := T\ t >oM t . The space M can be identified with the largest possible family of 
candidate test functions for weak solutions to either the stochastic heat equation ( 12. 31) or the 
stochastic wave equation (12.41) . Standard heuristics from PDEs then tell us that (12. 3p and/or 
(12.41) has random-field solutions if and only if 5 X G M for all x G R^; this can be interpreted 
as an equivalent definition of random- field solutions. When 5 X G M we may write u(t ,x) in 
place of u(t , 5 X ) . In order to prove Theorem 12.11 we will need some a priori estimates on the 
weak solutions of both the stochastic equations (12.31) and (12. 4p . We proceed by identifying 
M with generalized Sobolev spaces that arise in the potential theory of symmetric Levy 
processes. Now let us begin by studying the parabolic case. 



3.1. The parabolic case. 

Proposition 3.1. Let H denote the weak solution (12. 9p to the stochastic heat equation (12.31) . 
Then, for all <p G y(R d ), A > 0, and t > 0, 

1 _ p -2t/A P 2t/X 

(3.2) S{\ ;<p) < E (\H(t , V )\ 2 ) < — *(A ; <p), 

where 

<3 ' 3) '"WfL (l/A) +Re^(fl ^ 

Next we record the following immediate but useful corollary; it follows from Proposition 
13. II by simply setting A := t. 
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Corollary 3.2. If H denotes the weak solution to the stochastic heat equation (\2.3\i . then 
for all ip e y(R d ) and t > 0, 

(3.4) K(^) <E(\H(t, V )\ 2 ) <4<f(t;^), 

The preceding upper bound for E(\H(t , y?)| 2 ) is closely tied to an energy inequality for the 
weakly asymmetric exclusion process. See Lemma 3.1 of Bertini and Giacomin (1999); they 
ascribe that lemma to H.-T. Yau. 

The key step of the proof of Proposition 13. II is an elementary real- variable result which we 
prove next. 

Lemma 3.3. If g : R + — > R + is Borel measurable and nonincreasing, then for all t, A > 0, 

poo rt /*oo 

(3.5) (l-e- 2 '/ A ) / e- 2s/x g(s)ds< g(s)ds<e 2t/x e - 2s/x g(s) ds. 

Jo Jo Jo 

Monotonicity is not needed for the upper bound on f Q g(s) ds. 

Proof of Lemma \3. 3\ The upper bound on J Q * g(s) ds follows simply because e 2 ^ s ^ x > 1 
whenever t > s. In order to derive the lower bound we write 



poo °° p(n+l)t 

/ e- 2s / x g(s)ds = J2 e- 2s / x g{s)ds 

JO n=Q Jut 

oo „ t 

< V e- 2nt/x / g{s + nt) ds. 

„-n Jo 



(3.6) 



Because g is nonincreasing we can write g(s + nt) < g(s) to conclude the proof. □ 
Proof of Proposition \3.1\ We know from (I2.9p and the Wiener isometry that 



(3-7) E(\H(t^)\ 2 )= fllPMl^ds. 



Since the Fourier multiplier of P* is exp(— s\l/(— £)) at £ G R d , we can apply the Plancherel 
theorem and deduce the following formula: 

(3-8) \\PMl HR * } = j^y d J Rd e- 2 ^M0\ 2 dt. 

Because Re\I/(£) > 0, Lemma [3.31 readily proves the proposition. □ 

Equation (13.31) can be used to define (f(A ; p) for all Schwartz distributions <p, and not 
only those in ,5^(R d ). Moreover, it is possible to verify directly that (p \— > <^(A ; ip) 1 ^ 2 defines 
a norm on J?(R d ). But in all but uninteresting cases, ^(R d ) is not complete in this norm. 
Let Z/^,(R d ) denote the completion of J^(R d ) in the norm S"{X ; •) 1 ^ 2 . Thus the Hilbert space 
L 2 ^(R d ) can be identified with M. The following is a result about the potential theory of 
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symmetric Levy processes, but we present a self-contained proof that does not depend on 
that deep theory. 

Lemma 3.4. The space L^,(R d ) does not depend on the value of X. Moreover, L^(R d ) is 
a Hilbert space in norm $(\ ; •) 1 / 2 for each fixed A > 0. Finally, the quasi-isometry (13.21) is 
valid for all t > , A > ; and (p G L 2 ^ 

Proof. We write, temporarily, Z/^ A (R d ) for L\ (R d ), and seek to prove that it is independent 
of the choice of A. 

Define for all distributions ip and i/j, 



(3.9) 



2(2tt) 



R „ (l/A) + He*({) ; RJ (1/A) + Re*({) 



For each A > fixed, ((p,ip) i— > S{\ ;ip is a pre-Hilbertian inner product on <$ fi (R, d ), and 

Thanks to Proposition 13.11 for all a > there exists a finite and positive constant c = 
c aj A such that c~ x <§ [a ; ip) < ${\\ip) < c<S{a\p>) for all p> G J7(Tl d ). This proves that 
L^ x {H d ) = L^ a (R d ), whence follows the independence of L^(R d ) from the value of A. 
The remainder of the lemma is elementary. □ 

The space L^(R d ) is a generalized Sobolev space, and contains many classical spaces of 
Bessel potentials, as the following example shows. 

Example 3.5. Suppose J? = — (— A) s / 2 for s E (0,2]. Then, Jzf is the generator of an 
isotropic stable-s Levy process, and L%(R d ) is the space H_ s /2(R- d ) of Bessel potentials. For 
a nice pedagogic treatment see the book of Folland (119761 Chapter 6). □ 

3.2. The hyperbolic case. The main result of this section is the following quasi-isometry; 
it is the wave-equation analogue of Proposition 13.11 

Proposition 3.6. Suppose the symmetry condition (I2.12p holds, and let W := {W(t ,tp); t > 
0, <p G J^(R d )} denote the weak solution to the stochastic wave equation (12. 4p . Then, 

(3.10) \t£ {t 2 -p) <V{\W{t ,p)\ 2 ) <2t£ {t 2 -^) ■ 

for all t > and p> G =5^(R d ). Moreover, we can extend W by density so that the preceding 
display continues to remain valid when t > and p> G L^(R d ). 



Proof. Although Lemma 13731 is not applicable, we can proceed in a similar manner as we did 
when we proved the earlier quasi-isometry result for the heat equation (Proposition 13. ip . 
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Namely, we begin by observing that 



i rt r Sin 2 S I 

(3.1D B(iKr (t , v) |») = _y o y— > \m\^ 

See (12.211) . If 9 > then sin# is at most the minimum of one and 9. This leads to the 
bounds 

(312) ^w^<^[J^^)\m?«*. 

£ (4L( f2A 4) 1^''* 

The upper bound follows from this and the elementary inequality t 2 A z^ 1 < 2/(t~ 2 + z), 
valid for all z > 0. 

In order to derive the [slightly] harder lower bound we first rewrite (13.111) as follows: 

(3.13) E(mt, V )\ ) - ^pi RJ ^- 2Vmt J -*(|)- * 

We shall analyze the integral by splitting it according to whether or not \I/ < 1/t 2 . 

Taylor's expansion [with remainder] reveals that if 9 is nonnegative, then sin# is at most 
9 - (9 3 /Q) + (# 5 /120). This and a little algebra together show that 

sin^ 29 2 

(3.14) 1 — > if < 6> < 2. 

Consequently, 

For the remaining integral we use the elementary bound 1 — (sin 6/ 9) > 1/2, valid for all 
9 > 2. This leads to the following inequalities: 



(3.16) 



, / sin (2^W) t) \ m){ 2 i r {m 

*(£) 



2 

-de 



2 w,v tA ^) ,l ^ )|de - 



The proof concludes from summing up equations (13.151) and (13.161) . and then plugging the 
end result into (I3.13p . □ 
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We now give a proof of Theorem 12.11 

Proof. Let us begin with the proof in the case of the stochastic heat equation (12.31) . Propo- 
sition 13.11 is a quasi-isometry of the maximal space M of test functions for weak solutions 
of (I2.3p into Z/^,(R d ). Thus, M can be identified with the Hilbert space Z/^,(R d ), and hence 
(12. 3p has random-field solutions if and only if S x G Ll^>(R d ) for all x G R d . Thanks to Lemma 
13.41 the stochastic heat equation (I2.3P has random-field solutions if and only if 

(3.17) / - — — — 7-7 < oo for some, and hence all, $ > 0. 

The first part of the proof is concluded since condition ( 13. 17ft is known to be necessary as well 



as sufficient for X to have local times (Hawkes, 1986, Theorem 1). The remaining portion of 



the proof follows from the preceding in much the same way as the first portion was deduced 
from Proposition 13.11 □ 



4. Spatial continuity: Proof of Theorem 12.41 

Proof. We work with the stochastic heat equation (12.31) first. Without loss of generality, 
we may — and will — assume that (12. 3p has a random- field solution H(t,x), and X has local 
times. Else, Theorem 12.11 finishes the proof. 

Let <f :— 5 X — 5 y , and note that |y3(£)| 2 = 2(1 — cos(£(x — y))) is a function of x — y. 
Because H(t , <p) = H(t , x) — H(t , y), equations (13.71) and (I3.8P imply that z t— > H(t , z) is a 
centered Gaussian process with stationary increments for each fixed t > 0. 

Consider the function 

1 r°° 1 _ cos(r£) 

(4.1) h(r) := — / - ^ y dC, defined for all r > 0. 

2tt y_ 00 1 + Rew(4) 

According to Lemma E31 Proposition [3JJ holds for all Schwartz distributions if G L 2 ^ (R d ). 
The existence of random-field solutions is equivalent to the condition that 5 X G L^(R d ) for 
all x G R. We apply Proposition 13. II to ip := 5 X — 5 y , with A := 1 [say], and find that 

(4.2) (1 - e~ 2t ) h{\x-y\) < E (\H(t , x) - H(t , y) | 2 ) < e 2t h(\x -y\). 
Define h to be the Hardy-Littlewood nondecreasing rearrangement of h. That is, 

(4.3) h(r) := inf{y > : g(y) > r} where g(y) := meas {r > : h(r) < y} . 

Then according to the proof of Corollary 6.4.4 of Marcus and Rosen (2006, p. 274), the 
stationary-increments Gaussian process x t— > H(t , x) has a continuous modification iff 

(4-4) / -Snmdr <oo. 



o+ 



f | log T | 
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Next we claim that (I4.4p is equivalent to the continuity of the local times of the symmetrized 
Levy process X. We recall that the characteristic exponent of X is 2Re and hence by the 
Levy-Khintchine formula it can be written as 

/oo 
(1 -co8(£x))v(dx), 
-oo 

where v is a cr-finite Borel measure on R with (1 A x 2 ) u(dx) < oo. See, for example 
Bertoin (119961 Theorem 1, p. 13). 
Suppose, first, that 

/CO 
(1 A \x\) u(dx) = oo. 
-oo 

Then, (14.41) is also necessary and sufficient for the [joint] continuity of the local times of the 
symmetrized process X; confer with Barlow (1988, Theorems B and 1). On the other hand, 
if (14.61) fails to hold then X is a compound Poisson process. Because X is also a symmetric 
process, the Levy-Khintchine formula tells us that it has zero drift. That is, X t = Y^=i ^i, 
where are i.i.d. and symmetric, and II is an independent Poisson process. It follows 

immediately from this that the range of X is a.s. countable in that case. This proves that the 
occupation measure for X is a.s. singular with respect to Lebesgue measure, and therefore 
X cannot possess local times. But that contradicts the original assumption that X has local 
times. Consequently, ( 14.411 implies ( 14. 6p . and is equivalent to the spatial continuity of [a 
modification of] local times of X. This proves the theorem in the parabolic case. 

Let us assume further the symmetry condition ( 12. 12ft . Because S'it 2 ] ip)/S{l ; <p) is bounded 
above and below by positive finite constants that depend only on t > (Proposition 13. ip . we 
can apply the very same argument to the stochastic wave equation (12.41) . but use Proposition 
13.61 in place of Proposition 13.11 This completes our proof of Theorem 12.41 □ 

5. An aside on temporal continuity 

We spend a few pages discussing matters of temporal continuity — especially temporal 
Holder continuity — of weak solutions of the stochastic heat equation ( 12.31) . as well as the 
stochastic wave equation (12.41) . 

Definition 5.1. We call a function g(s) a gauge function if the following are satisfied: 

(1) g(s) is an increasing function; 

(2) g(s) is a slowly varying function at infinity; 

(3) g(s) satisfies the integrability condition, 

(5.1) / - - f. - . . < oo. 



o+ 



s\og(l/s)g(l/s) 
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Next, we quote a useful property of slow varying functions (Bingham et al., 1989, p. 27). 

Proposition 5.2. If g is a slowly varying function and a > 1, then the integral t~ a g{t) alt 
converges for every x > 0, and 

5.2 / »±±dt~- — T as x^oo. 

K J J x t a (a-l)a— 1 

We can now state the main theorem of this section. It gives a criteria for the temporal 
continuity of the weak solutions of our stochastic equations. 

Theorem 5.3. Let H denote the weak solution to the stochastic heat equation (12.31) . Let 

g be a gauge function in the sense of Definition 15. 1\ Choose and fix tp G L"y(R, d ). Then, 
1 1 — > H(t , ip) has a continuous modification if the following is satisfied 

r*i\ f Mi + Refr(QMi + |vi/(0|) ,^, 2 ^ 

(5 ' 3) L 1 + Re*(0 di < °°- 

Moreover, the critical Holder exponent o/t i — > H(t,(p) is precisely 

(5.4) md*(.;<p) := hminf l ° g f { " ; v) . 

40 loge 

Consequently, t \— > H(t,<p) has a Holder- continuous modification (a.s.) iff ind ${• ; (p) > 0. 

// the symmetry condition (12.121) holds, then f !5.3j) guarantees the existence of a continuous 
modification of t i— > W(t,</?) ; where W denotes the weak solution to the stochastic wave 
equation. Furthermore, (15 .4p implies the temporal Holder continuity oft \— > W(t , <£>) of any 
order < | ind ; </?) . 

We now give two examples. The first one is about the temporal Holder exponent while the 
second one provides a family of random-field solutions which are almost surely unbounded 
in every open space-time set. 

Example 5.4. Suppose d = 1 and Jzf = —(—A)"/ 2 for some a G (1,2]. Choose and fix 
x G R. According to Theorem 12 .![ 5 a G //^(R), so we can apply Theorem 15.31 with ip := 5 X . 
In this case, 

1 f°° dC, 

(5.5) ^ (eA ' ) = ^/-oo(lA) + ie 

= const -e 1 "^. 

In particular, ind &{• ;8 X ) = 1 — (1/a), whence it follows that the critical Holder exponent of 
t i— > H(t , x) is precisely | — (2a) -1 . When j£f = A, we have a = 2 and the critical temporal 
exponent is 1/4, which agrees with a well-known folklore theorem. For example, Corollary 
3.4 of Walsh (I1986[ pp. 318-320) and its proof contain this statement for the closely-related 
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stochastic cable equation. In particular, see the last two lines on page 319 of Walsh's lectures 
(loc. cit.). □ 



Example 5.5. Choose and fix a G R. According to the Levy-Khintchine formula ( |Bertoin, 1996 
Theorem 1, p. 13) we can find a symmetric Levy process whose characteristic exponents sat- 
isfies ^(O ~ £(log£) a as £ — > oo. Throughout, we assume that a > 1. This ensures that 
condition fl3TTI) is in place; i.e., (1 + tf)" 1 G L X (R). Equivalently, that both SPDEs (TO) 
and (12.41) have a random-field solution. 

A few lines of computations show that for all x G R, 



(5.6) d?(e;<y x|log(l/ 



-a+l 



where a(e) x 6(e) means that a(e)/6(e) is bounded and below by absolute constants, uni- 
formly for all e > sufficiently small and x G R. Next, consider 

(5.7) d(s , t) := yfE(\H(s,x)-H(t,x)\ 2 ). 
According to Proposition 15.61 below, (15. 6ft implies that 

/ . \ (l-a)/2 

(5.8) logf-— 

uniformly for all s and t in a fixed compact subset [0,T] of R + , say. Let iV denote the 



metric entropy of [0 , T] in the [pseudo-] metric d ( Dudley, 19 67). That is, for all e > 0, we 
define N(e) to be the minimum number of d-balls of radius e needed to cover [0 , T]. Then, 
it is easy to deduce from the previous display that log N(e) x e~ 2 ^ a ~ l \ and hence 

(5.9) lime v / logiV(e) = if and only if a > 2. 

40 



Therefore, if 1 < a < 2 then Sudakov minorization (Marcus and Rosen, 2006, p. 250) tells 
us that the process t i— > H(t,x) — and also t h- > W(t,x) — does not have any continuous 
modifications, almost surely. On the other hand, if a > 2, then the integrability condition 
(15. 3ft holds manifestly, and hence t i— > i/(t , x) and t i— > (t , x) both have continuous 
modifications. Thus, the sufficiency condition of Theorem [531 is also necessary for the present 
example. In addition, when a < 2, the random-field solutions H and W are both unbounded 
a.s. in every open set. This assertion follows from general facts about Gaussian processes; 
see, for example equation (6.32) of Marcus and Rosen ( 2006 , p. 250). We skip the details. □ 

5.1. Estimate in the parabolic case. 

Proposition 5.6. For every t,e>0 and (p G L^(R d ), 

(5.10) \ S{e ■ ip) < E (\H(t + e,(p)-H(t, V )\ 2 ) < <f (e ; <p) + e 2t ^(e ; <p), 
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where 



("1) , {c ., vh= ^f 



Proof. Because of density it suffices to prove that the proposition holds for all functions 
(p G J^(R d ) of rapid decrease. By (12.91) and Wiener's isometry, 

E(\H(t + e^)-H(t, V )\ 2 ) = f ( \(P t *_ s+e <p)(y)-(P t *_ sV )(y)\ 2 dyds 
(5.12) Jo jRd 



t+e 
t jR d 



(Pt-s+eV) (y)\ dyds. 



We apply Plancherel's theorem to find that 

E(\H(t + e,<p)- H{t , V )\ 2 ) = JL jf ' jf ^ |e- (s+e) * ( -« - e~^\ 2 \0(O\ 2 ds 



(5.13) 



1 



(27r) d 

Thus, we can write 



e-^^\m\ 2 dyds. 

^R d 



(5.14) E(\H(t + e,<p)-H(t,<p)?) : = 
where 

(5.15) T i: =r/ e- 2sRe *«)|l-e- £ ^| 2 |^(0r^^, 

Jo Jn d 

and 

(5.16) T 2 := f / e- 2sRe ^|^(0| 2 ^^. 
First we estimate T 2 , viz., 

(5.17) / e~ 2sRe ^ds = e — -—. 

1 1 Jo 2eRe*(0 

Because 

11 1 — p -61 2 

(5.18) < — < for all 6 > 0, 

v; 21+fl-fl-l+fl 

it follows that \(2ir) d g{e\y) <T 2 < (27c) d £{e; <p). Since T\ > we obtain the first inequality 
of the proposition. 

According to Lemma [3.31 with A = 1, 

e2t f h _- e ff(oi 2 l£(OI 5 



(5.19) T x <— i- e - eV W v 
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Because |1 - e - £ * (€) | 2 < 1 Ae 2 |^(£)| 2 , it follows that T\ < {27r) d exp(2t)J?(e ; <p), and hence 
the proof is completed. □ 



5.2. Estimate in the hyperbolic case. 



Proposition 5.7. Assume the symmetry condition (j2.12p . and let W denote the weak solu- 
tion to the stochastic wave equation (12.41) . Then, for all t > , e > 0, and tp> e L 2 ^ (R d ), 



(5.20) 



E(\W{t + e,cp) -W(t,tp)\ 2 ) < (8t + 6e)^(e 2 ;^). 



Proof. By density, it suffices to prove the proposition for all functions ip G S^(R d ) of rapid 
decrease. Henceforth, we choose and fix such a function (p. 
In accord with (12.211) we write 

E (\W(t + e,tp) — W(t , <p)\ 2 ) := 7\ + T 2 , 



(5.21) 
where 

(5.22) 
and 

(5.23) 



(2vr) d 



T> 



sin(^(0(r + e)) - sin ■ 



^R d 



1 



*(£) 

€ /• sin 2 (a/^TO r 



(2vr) d 7 t 7 Rd *(£) 



|fp(e)| 2 dedr. 



We estimate T 2 first: The argument that led to (13 . 1 6[) also leads to the following inequality: 
sin 2 (V#(£)r) r f x 

R d 



(5.24) 



R rf 



Because J 2e (r 2 A a) dr < 3e (e 2 A a) for all a, e > 0, 



(5.25) 



To < 



3e 



e 2 A 



*(0 



1^(0 1 2 de 



< 6e<?(e 2 ;^). 

The estimate for Ti is even simpler to derive: Because | sin a — sin/?| 2 is bounded above 
by the minimum of 4 and 2[1 — cos(/3 — a)], 

sin (v^(0 (r + e)) - sin : 



(5.26) 



R" 



*(0 



< 



MOI 2 ^ 

1 - cos (v^(0 



Rrf 



*(£) 
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This and the elementary inequality 1 — cosx < x 2 /2 together yield the bound 



l-cos(v^(O e 



1 



(5.27) ±j- ^<2(e 2 A T ,. 

Consequently, 



1 



sin (y^(£) (r + e)) - sin (v^IO r 



2 



(5.28) p^y RJ v - ife ^ M^K)| 2 <«<8^), 

whence T\ is at most t times the right-hand side of the preceding. This and (I5.25P together 
yield the proof. □ 

Proof of Theorem \5.3l We start with the weak solution H to the stochastic heat equation. 
Throughout, <p G L\ (R d ) is held fixed. 

If the integrability condition (15.31) holds, then according to Lemma [5.61 for all T > 0, 

,, 0Q , E(\H(t + e, ¥ >)-H(t, ¥ >)\ 2 ) 

(5.29) sup — — - < 4e + 1 < oo. 

t,e>o-. & (e ; if) 

0<t<t+e<T 

Since e i— > ${e • ip) is nondecreasing, a direct application of Gaussian-process theory implies 
that {H(t ,<p)}te[o,T] has a continuous modification provided that 

(5.30) / V^2Z&<oo. 



See Lemma 6.4.6 of Marcus and Rosen (2006, p. 275). A standard measure-theoretic ar- 
gument then applies to prove that t i— > H(t , tp) has a continuous modification. A similar 
argument works for the weak solution W to the stochastic wave equation 02.4)) . but we appeal 
to Proposition 15.61 in place of 15.71 Thus, the first portion of our proof will be completed, 
once we prove that condition (15.31) implies (15.301) . 
Let us write 



(5.31) io+ ev /log(l/e) 7 + ev /log(l/e) 7 + e^/log(l/e) 

:= h + h- 

Let us consider l x first. We multiply and divide the integrand of I\ by the square root of 
g(l/e), and then and apply the Cauchy-Schwarz inequality to obtain the following: 



(5.32) 



o+ e y V.io+ elog(l/e)^(l/e) 

const . ( / £(^MV£), f y /2 . 

o+ e J 
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Note that 



I * = M ®\ « ^ 1 + eRe9(0 * 

/ 3 + /4 

' (2vr) d ' 

with the notation being clear enough. We will look at J 3 and I4 separately. We begin with 
I4 first. 

•W<Jo l + eRe^(0 
(5,34) < / \m\ 2 d£' [ 1/e 9(l/e)de 



JRc^<e JO 

< const ■ <§{\ ; (p), 

where we have used Proposition 15.21 to obtain the last inequality. In order to bound J 3 we 
let N > e and look at the following: 

ff(l/e) 



1 + eiV 



de 



(5.35) _ r /K ioz£i d£+ p-mm.* 



o 1 + eiV J 1/iV 1 + eiV 
= h + h- 

On one hand, some calculus shows that 

(5.36) 7. < «>. 

On the other hand, we can change variables and appeal to Proposition 15.21 and deduce that 



(5.37) 



Jn u 



N 

We obtain the following upon setting N := Re and combining inequalities fl5.34l) - fl5.36l) 

above: 
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We now look at l^- Let us recall the definition of ; <p) and write 
J?(e;y?) 1 



(■•■)*+/ (•■•)*+/ (•■■)« 

|*|<l/e ■'l/e<|4 r |<l/e J|*|>l/e 



(5.39) ' ' " ' ( 27r ) d 

= const ■ [J 7 + 7 8 + 2g] • 

Because J 7 < e 2 c?(i ; </?) whenever e < 1/e, 

(5.40) / 7 ;== de < const • ^(1 ;¥ >). 

Jo e v /log(l/e) 

An application of Cauchy-Schwarz inequality yields 



o ev /log(l/e) 7o VA/ e <|*|<iA 1 + Re*(0 J y/hg(l/e) 

(5 ' 41) <»■*■ (7 

\JkJo l + Re*(0 log(l/e) 

< const • a/ <§{\ ; y?). 

We multiply and divide the integrand below by the square root of g(l/e) and apply the 
Cauchy-Schwarz inequality in order to obtain 

1/2 



(5.42) 



o ev /log(l/e) y|*|>i/ e e(l + Re*(0) 



1/2 



In the last inequality we have changed the order of integration and used Proposition 15.21 
Taking into account inequalites (I5.36l) -( l5.42l) . we obtain 

Inequalities (15.381) and (15.431) together with (15.311) imply that 
(5.44) / vg5^SM d£ <oc, 



Jo+ ev/log(l/e) 

provided that ( 15. 3ft holds. This proves the first assertion of the theorem. 
Suppose 7 > 0, where 

(5.45) 7 :=ind<f (•;</?), 
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for brevity. Then by definition, £(e',<p) < e 7+ °^ 1 - ) as e | 0. Another standard result from 
Gaussian analysis, used in conjunction with Proposition 15.61 proves that H has a Holder- 



continuous modification with Holder exponent < 7/2 ( Marcus and Rosen, 2006 , Theorem 
7.2.1, p. 298). The proof for W is analogous. 

For the remainder of the proof we consider only the weak solution H to the stochastic heat 
equation ( 12. 3ft . and write H t := H(t,ip) for typographical ease. If 7 > 0, then, Proposition 
15.61 and elementary properties of normal laws together imply that 

(5.46) inf \\H t+e - H t \\ L2(p) > e 7+o(1) for infinitely-many e j 0. 



Consequently, for all 5 G (0 , 1), q > 7/2, T > S > 0, and t G [S , T] — all fixed — the following 
holds for infinitely- many values of e j 0: 



(5.47) 



p\ sup \H r+e -H r \ >e q \ >p\\H t+e -H t \ > 6 \\H t+e - H t \\ L2(P) \ 

[re[S,T] J L J 

= 1 - ^ J exp(-x 2 /2) dx. 
The inequality exp(— x 2 /2) < 1 then implies that 

(5.48) P \ sup \H r+e - H r \ > (1 + o(l))e q for infinitely many e J, \ > 1 - J^- 5. 

[r G [S,T] 

Since g > 7/2 is arbitrary, we can enlarge it arbitrarily to infer that 



\H r+t — H r \ I . _ 12 



(5.49) P < limsup sup = oo>>1 — \ — 5. 

\ 40 r6[S>] I V 7T 

Let 5 J, to find that 

/ \ i I ^v-\-e — Hf I 

(5.50) limsup sup = 00 a.s. 

40 re[S,T] £ q 

This proves that any q > 7/2 is an almost-sure lower bound for the critical Holder exponent 
of H. Moreover, in the case that 7 = 0, we find that (I5.50p holds a.s. for all q > 0. Thus, 
it follows that with probability one, H has no Holder-continuous modification in that case. 
The proof is now complete. □ 



6. Spatial and joint continuity: Proofs of Theorems 12.51 and 12.61 

Proof of Theorem \2.5[ We begin by proving the portion of Theorem 12.51 that relates to the 
stochastic heat equation and its random-field solution {H(t , x); t > 0, x G R}; Theorem 12. II 
guarantees the existence of the latter process. 
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Throughout we can — and will — assume without loss of generality that x i— ► H(t , x) is 
continuous, and hence so is the local time of X in its spatial variable. As we saw, during the 
course of the proof of Theorem 12.41 this automatically implies the condition ( 14. 6ft . which we 
are free to assume henceforth. 



According to Proposition 13.11 specifically Corollary | 

(6.1) \£(i ; 8 X - 8 y ) < E (\H(t , ar) - H(t , y) | 2 ) < 4<£(t ; 6 X - 5 y ). 
Since |4(£) - 4(0 1 2 = 2 [ 1 _ cos(£(x - y))], it follows from this and (H~Tj) that 

(6.2) <E(|tf(t,;c)-tf(i,y)| 2 ) <8/i(|z-j/|). 

This implies that the critical Holder exponent of z i— ► H(t ,z) is almost surely equal to 
one-half of the following quantity: 

log/i(e) 



(6.3) ind h := liminf . 

e|o loge 

We will not prove this here, since it is very similar to the proof of temporal Holder continuity 
(Theorem 15.31) . Consequently, 

(6.4) z t— > H(t , z) has a Holder-continuous modification iff ind h > 0. 

Among other things, this implies the equivalence of parts (1) and (2) of the theorem. 

Let Z(t,x) denote the local time of X at spatial value x at time t > 0. We prove that 
(1), (2), and (3) are equivalent by proving that (16.41) continues to hold when H is replaced 
by Z. Fortunately, this can be read off the work of Barlow (I1988p . We explain the details 
briefly. Because ( 14.61) holds, Theorem 5.3 of Barlow 019881) implies that there exists a finite 
constant c > such that for allt > and finite intervals I C R, 

ffitf r \Z(t,a)-Z(t,b)\ ( A 1/2 

6.5 hm sup — > c sup Z(s , x) a.s. 

<U° a,bei v/M|a-6|)log(l/|6-a|) ~ \x/i J 

|a-b|<<5 

Moreover, we can choose the null set to be independent of all intervals I C R with rational 
endpoints. In fact, we can replace / by R, since a \— > Z(t , a) is supported by the closure 
of the range of the process X up to time t, and the latter range is a.s. bounded since X is 
cadlag. 

By their very definition local times satisfy Z(s ,x)dx = s a.s. Thus, sup^gR Z(t , x) > 
a.s., whence it follows that for all q > \ ind h, 

(r. p\ r \Z(t,a)-Z(t,b)\ 

(6.6) hm sup : = cxo a.s. 

*i° \a-b\<6 \a - b\i 
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That is, there is no Holder-continuous modification of a i— ► Z(t,a) of order > | ind h. In 
particular, if ind h = 0, then a i— > Z(t , a) does not have a Holder-continuous modification. 

Define d(a , 6) := a/MI° — fej); it is easy to see that d is a pseudo-metric on R. According 
to Bass and Khoshnevisan (I1992p . 

(R7\ v \Z(t,a)-Z(t,b)\ ( A 1/2 

(0.7) iimsup sup ,, < 2 supZ t.s a.s., 

«5i0 d(a, b )<6 f d(a > b) (log N(u) y/*du ~ V/ 7 

where N(u) denotes the smallest number of <i-balls of radius < u needed to cover [—1 , 1]. This 
sharpened an earlier result of Barlow (1985, Theorem 1.1). Furthermore, sup x Z(t ,x) < oo 
a.s. (B ass and Khoshnevisan, 1992 Theorem 3.1). Consequently, 

(6.8) sup \Z(t,a)- Z(t,b)\=o( [ (log N(u)) 1/2 du] as 5 | 0. 

d{a,b)<5 \Jo J 

Since h is increasing, 

(6.9) sup \Z(t,a)-Z(t,b)\ = 0\ (log N(u)) 1/2 du) as 5 | 0. 

|o-6|<<5 \Jo J 

According to equations (6.128) and (6.130) of Marcus and Rosen ( 12006^ Lemma 6.4.1, p. 
271), there exists a finite constant c > such that for all u > small, 

N(u) < 



£ 2 {(x,y)G[-l,l] 2 : h(\x-y\)<u/A} 
^ ' ' const 

- h-\ u /Ay 

where ^ 2 (^4) denotes the Lebesgue measure of A C R 2 . [Specifically, we apply Lemma 6.4.1 
of that reference with their K := [—1 , 1] and their /i 4 := c] This and (16. 9p together imply 
that with probability one the following is valid: As 5 j 0, 

sup \Z(t,a)-Z(t,b)\=o( f | log«| 1/2 /*((///) 
(6.11) la - bl<5 

= n (hfxvirxrh /^l 1 /2^ _u n ( 

l wllogul 1 / 2 

The last line follows from integration by parts. If 7 := ind h > 0, then h(5) = o(5 q ) for all 
fixed choices of q G (0,7/2). It follows that a t— > Z(t,a) is Holder continuous of any order 
< 7/2. Among other things, this implies (\6A\1 with Z replacing H, whence it follows that 
(l)-(3) of the theorem are equivalent. 

The hyperbolic portion of the theorem is proved similarly, but we use Proposition 13.61 in 
place of Proposition 13.11 everywhere. □ 



LOCAL TIMES AND SPDES 



27 



Proof of Theorem \2.b\ Since f3" is positive, Re\l/(£) — > oo as |£| — > oo. Therefore, for all 
e, $ > and x G R, 

^ ^ x)= (2^F L (l/e)+1te*(0 
< const • e - 



(l/e)+Re^(0' 

If £ is sufficiently large, then for all 7 G (d,P"), we can find a constant C 7 G (0 , 00) such 
that Retf (0 > C 7 |£| 7 for all £ G {Re^ > ■&}. Consequently, 

(6.13) V 7 { Rc*>.}(lA) + |elV 

= O (e 1 "^) , 

as e I 0. Thus, ind<^(» ; 8 X ) < 1 — (d/(3"), where this index was introduced in ( 15.4h . A similar 
calculation shows that ind h < (3" — d; confer with (16.31) for the definition of this quantity. 
Thus, for all fixed T > 0, Proposition 15.61 and (16.21) together prove the following: For all 

r < (p" - d)/P", ( < P" - d, x, y G R, and s, t G [0 , T], 

(6.14) E(\H(t,x) - H(s,y)\ 2 ) < const ■ (\s - t\ T + \x - y\ c ) . 

A two-dimensional version of Kolmogorov's continuity theorem finishes the proof; see the 
proof of Theorem (2.1) of Revuz and Yor (1991, p. 25). 

The proof, in the case of the stochastic wave equation, is similar, but we use Propositions 
13.61 and 15.61 instead of Propositions 13.11 and 15.71 respectively. □ 

7. Heat equation via generators of Markov processes 

We now consider briefly the stochastic heat equation, where the spatial movement is gov- 
erned by the generator ££ of a [weakly] Markov process X := {X t }t>o that takes values in a 
locally compact separable metric space F. We assume further that X admits a symmetrizing 
measure m that is Radon and positive. Let us emphasize that m satisfies (Ptf, g) = (/, Ptg) 
for all t > and /, g G L 2 (m), where {Pt}t>o denotes the transition operators of X, and 
((fi , (f 2 ) := J (f\(p2 dm for all ip 2 G L 2 (m). 

7.1. The general problem. Consider the stochastic heat equation 

d t u(t , x) = {££u){t , x) + w(t , x), 
u(0,x) = 0, 

valid for alH > and x G F . Here, the underlying noise w in (17. ip is a Gaussian martingale 
measure on R + x F in the sense of Walsh ( ]1986j) : w is defined on the filtered probability 
space (f2 , J^~, {^t}t>o, P) an d w t {tp) '■— f J F (f(s , x) w{dx ds) defines an {J^} 4 > -martingale 



(7.1) 
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for ip G L 2 (ds x m); w can be characterized by the covariance functional for the corresponding 
Wiener integrals: 

(7.2) E (^J fdw ■ J gdw^j = J J J f(s,x)g(s,y)rn(dx)m(dy)ds, 

for all /, g G L 2 (ds x m). 

We can follow the description of Walsh (1986) and write the weak form of equation (17. II' 
as follows: For all ip G L 2 (m) and £ > 0, 

f 



(7.3) u(t,(p) = / / (P t - S <p)(x)w(dxds). 

JO J F 

Lemma 7.1. The integral defined by (17.31) is well defined for all ip G L 2 (m) 



Proof. We follow closely the proof of Proposition 12.81 and apply the fact that the semigroup 
{P s } s > is a contraction on L 2 (m). □ 

Let Z denote the occupation measure of X; consult (11.31) . The following is the key result 
of this section. It identifies an abstract Hilbertian quasi-isometry between the occupation- 
measure L 2 -norm of X and a similar norm for the solution to the stochastic heat equation 
(JUD for J5f. 

Theorem 7.2. If u denotes the weak solution to (17. ip . then for all p G L 2 (m) and t > 0, 

(7.4) | tE (\u(t , p)\ 2 ) < E m (\Z(t , ^)| 2 ) < 4tE (\u(t , p)\ 2 ) . 

As usual, E m refers to the expectation operator for the process X, started according to 
the measure m. 

The preceding theorem follows from the next formula. 

Proposition 7.3. If u denotes the weak solution to (17. ip . then for all p G L 2 (m) and t > 0, 

(7.5) E m (\Z(t,p)\ 2 ) =4 fE(\u{s/2^)\ 2 ) ds. 

Jo 



Proof of Proposition 7.3. Since m is a symmetrizing measure for X, the P m -law of X u is m 



for all u > 0. By the Markov property and Tonelli's theorem, 



(7.6) 



E m (\Z(t 1 p)\ 2 ) = 2 / / (P v ^ u (p , p) dv du. 

Jo Ju 
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We computing the Laplace transform of both sides, viz., 

»oo pt ft 



poo poo pt pt 

/ e- xt E m (\Z(t,<p)\ 2 ) dt = 2 / / e~ xt (P v ^ u <p,<p) dvdudt 

JO JO Jo Ju 



oo /*oo 



(7.7) =2/ / (/ e-^tft ) (P v ^(p,(p) dvdu 

Jo Ju 

2 



oo 

As 



e- As (P s if , yj) (is. 



The exchange of the integrals is justified because [P r ^p , <p) = \\P r /2 ^ll^rm) * s positive and 
finite. Because 93 G L 2 (m), Fubini's theorem implies that for all A > 0, 



(7.8) 



e- M E m (\Z(t,if)\ 2 ) dt=-{Rx<P,<P) 



where R\ := / °° exp(— Xs)P s ds defines the resolvent of {-Pt}t>o- 

Let T\ denote an independent mean- (1/ A) exponential holding time. The preceding display 
can be rewritten as follows: 

(7.9) E m (|Z(T A ,<^)| 2 ) = ~{R x <p,<p) for all A > 0. 

Next we consider the weak solution u to the stochastic heat equation (17.11) by first observing 
that E(Kt,^)| 2 ) = j l Q \\PM\ i 2 (m)^ s - ^ follows from this, and successive applications of 
Tonelli's theorem, that for all (3 > 0, 

poo pt 

E(\u{T p ,<p)\ 2 )= fie-* I \\PM\h (m )dsdt 
Jo Jo 



(7.10) 



e- ps \\P s p\\ 2 L2{m) ds. 



Because H-fsVllz^fm) = (^s^jf), we may apply Fubini's theorem once more, and select 
p := 2A, to find that 

poo 

E(KT 2A ,<^)| 2 ) = / e- 2Xs (P 2s cf,if) ds 
Jo 

an) =(<p, r e - 2xs p 2s(P ds 



The condition of square integrability for <p justifies the appeal to Fubini's theorem. We can 
compare (17. 9p and (17.111) to find that 

(7.12) E m (\Z(T x ,<p)\ 2 ) = jE (\u(T 2X ,<p)\ 2 ) for all A > 0. 
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Define q(t) := E m (\Z(t , <^)| 2 ), p(t) := E(\u(t/2 , ^)| 2 ) and l(t) := 1 for all t > 0. The 
preceding shows that the Laplace transform of q is equal to 4 times the product of the 
respective Laplace transforms of p and 1. Thus, we can invert to find that q = 4p* 1, which 
is another way to state the theorem. □ 



Proof of Theorem \7.2\ Let us choose and fix a measurable function : F — > R such that 
| </?| G L 2 (m) . The defining isometry for Wiener integrals yields the following identity, where 
both sides are convergent: E(\u(t , y?)| 2 ) = / ||P s <^|| 2 2 / m ) <is. Proposition 17.31 implies that 

(7.13) 2tE(Kt/4,^)| 2 ) <E m (\Z(t, V )\ 2 ) <4tE(|n(t/2,^)| 2 ). 

[For the lower bound, we use the bound £E(\u(s/2 , <p)\ 2 ) ds > J t/2 E{\u{s/2 , y?)| 2 ) ds.} By 
monotonicity, E(\u(t/2 ,ip)\ 2 ) < E(\u(t , <f)\ 2 ), whence follows the announced upper bound 
for E m (\Z(t,<p)\ 2 ). 

In order to prove the other bound we first write 

(7.14) Z(t , if) = Z(t/2 , ip) + Z(t/2 , if) o 9 t/2 , 

where {6 s } s >o denotes the collection of all shifts on the paths of X. We care only about 
distributional properties. Therefore, by working on an appropriate probability space, we can 
always insure that these shifts can be constructed; see Blumenthal and Getoor fll968p . 
We apply the Markov property at time t/2. Since P m o X~, 2 = m, it follows that 

(7.15) E m (\Z(t/2 , ip) o 9 t/2 \ 2 ) = E m (\Z(t/2 , V )\ 2 o 6 t/2 ) = E m (\Z(t/2 , <^)| 2 ) , 
and hence E m (\Z(t , ip)\ 2 ) < AE m (\Z(t/2 , ip)\ 2 ). Consequently, 

(7.16) E m (\Z(t,v)\ 2 ) < 16E m (|Z(t/4,v?)| 2 ) for all t > 0. 

This and the first inequality of (IT. 13j) together imply the remaining bound in the statement 
of the theorem. □ 



7.2. The stochastic heat equation in dimension 2 — e. We now specialize the setup of 
the preceding subsection to produce an interesting family of examples: We suppose that F 
is a locally compact subset of R d for some integer d > 1, and m is a positive Radon measure 
on F, as before. Let {-Ra}a>o denote the resolvent of X, and suppose that X has jointly 
continuous and uniformly bounded resolvent densities {rA}A>o- m particular, r\(x,y) > 
for all x,y G F and 

(7.17) (Rxf)(x)= [ r x (x,y)f(y)m(dy), 

Jf 
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for all measurable functions f : F —> R + . Recall that X has local times {Z(t , x)} t >o, x &F if 
and only if for all measurable functions f : F —* R+, and every t > 0, 

(7.18) Z(t,f)= [ Z(t,z)f(z)m(dz), 

J F 

valid P x -a.s. for all x £ F . Choose and fix some point a £ F, and define 

(7.19) f?(z) := lB{a ; e) ^\, for all z £ F and e > 0. 

Of course, _B(a , e) denotes the ball of radius e about a, measured in the natural metric of 
F. Because r\ is jointly continuous, lim e i (R\f®)(x) = r\(x,a), uniformly for x-compacta. 
Define (p e ^ := / e a — fg, and observe that ip £ L 2 (m). Furthermore, 

(7.20) lim (R x cp e , 5 , <p e>S ) = ^A(Rxf: , ft) ~ 2 (Rxf? , / 5 a ) + {Rxff , / 5 a )} = 0. 

If /i £ L 2 (m), then the weak solution h to (17. 11) — where Jzf denotes the L 2 -generator of 
X — satisfies 

-t 



E(\u(t,h)\ 2 ) = [ \\P s h\\ 2 LHm) ds 
Jo 



(7-21) <e^ e-^\\P s h\\U m) ds 



•X- 



e 



L 2 (m) 

JO 

2Xt 



-(R x h,h). 



2 

Therefore, Theorem 17.21 and (17.201) together imply that {u(t , f^)} t >o is a Cauchy sequence 
in L 2 (P) for all t > 0. In other words, we have shown that the stochastic heat equation (17.11) 
has a "random- field solution." 



Example 7.4. It is now easy to check, using the heat-kernel estimates of Barlow (1998, 
Theorems 8.1.5 and 8.1.6), that for all d £ (0,2) there exists a compact "fractal" F C R 2 
of Hausdorff dimension d such that Brownian motion on F satisfies the bounded/continuous 
resolvent-density properties here. The preceding proves that if we replace Jzf by the Laplacian 
on F, then the stochastic heat equation (17.11) has a "random-field solution." Specifically, the 
latter means that for all t > and a £ F, u(t , ff) converges in £ 2 (P) as e [ 0, where /" 
is defined in (17.191) . This example comes about, because the fractional diffusions of Barlow 
(I1998p have local times when the dimension d of the fractal on which they live satisfies d < 2. 
See Barlow (119981 Theorem 3.32), for instance. □ 



U) 
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8. A semilinear parabolic problem 

We consider the semilinear problems that correspond to the stochastic heat equation (12.31) . 
At this point in the development of SPDEs, we can make general sense of nonlinear stochastic 
PDEs only when the linearized SPDE is sensible. Thus, we assume henceforth that d — 1. 

We investigate the semilinear stochastic heat equation. Let b : R —>■ R be a measurable 
function, and consider the solution H b to the following SPDE: 

d t H b (t,x) = (3?H b )(t,x) + b(H b (t,x))+w{t,x), 

H b (0 ,x) = 0, 

where =Sf denotes the generator of the Levy process X, as before. 

Equation (18.11) has a chance of making sense only if the linearized problem (12.31) has a 
random-field solution H, in which case we follow Walsh (1986]) and write the solution H b as 
the solution to the following: 

i*t roc 

(8.2) H b (t ,x) = H(t,x)+ / b(H b (s ,x-y)) P^ s (dy) da, 

where the measures {Pt}t>o are determined from the semigroup of X by Pt(E) := (P(ls)(0) 
for all Borel sets E C R. [This is standard notation.] We will soon see that this random 
integral equation has a "good solution" H b under more or less standard conditions on the 
function b. But first, let us make an observation. 

Lemma 8.1. // (12.31) has a random-field solution, then the process X has a jointly measur- 
able transition density {pt{x)}t>o,xen that satisfies the following: For all r\ > there exists 
a constant C := C v G (0 , oo) such that for all t > 0, 



(8.3) / \\ Ps \\ 2 L2{R) ds<Ce*. 

Jo 

Finally, (t , x) i— > Pt{x) is uniformly continuous on [e , T] x R for all fixed e, T > 0. 



Proof. We can inspect the function y = xexp(— x) to find that exp(— x) < (1 + x) 1 for all 
x > 0. Consequently, 

The second integral, however, has been shown to be equivalent to the existence of random- 
field solutions to (12.31) : see (13. 17ft . It follows that the first integral in (18 Ah is convergent. We 
apply the inversion theorem to deduce from this that the transition densities of X are given 
by Pt(x) = (27r) _1 exp{— ix^ — t\E r (^)}<i^, where the integral is absolutely convergent for 
all t > and x G R. Among other things, this formula implies the uniform continuity of 
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Pt(x) away from t = 0. In addition, by Plancherel's theorem, for all s > 0, 

1 /*oo 



oc 



( 8 - 5 ) i r - 

27T loo 

Therefore, Lemma [3.31 and Tonelli's theorem together imply that for all A > 0, 

-t p 2i/A /-oo ^ 



(8.6) y o wi»cR) ^ < 



< oo. 



(1/A)+Re^(0 

This completes the proof. □ 
Thanks to the preceding lemma, by (18. ip we mean a solution to the following: 

l*t i*oo 

(8.7) H b (t,x)=H(t,x)+ / b(H b (s,x-y))p t - s (y)dyds. 

JO J -oo 

For the following we assume that the underlying probability space (Q , J^", P) is complete. 

Theorem 8.2. Suppose b is bounded and globally Lipschitz, and the stochastic heat equa- 
tion (12.31) has a random-field solution H; thus, in particular, d = 1. Then, there exists a 
modification of H , denoted still by H , and a process H b with the following properties: 

(1) H b G L p loc (R + x R) for all p G [1 , oo). 

(2) With probability one, ( \8.7b holds for all (t , x) G R+ x R. 

(3) For all T > 0, J is a.s. bounded and continuous on [0 ,T] x R, where 

l*t /'OO 

J(t,x):= / b(H b (s ,x - y))pt- s (y) dyds. 

JO J -oo 

Remark 8.3. Before we proceed with a proof, we make two remarks: 

(1) It is possible to adapt the argument of Nualart and Pardoux (119941 Proposition 
1.6) to deduce that the laws of H and H b are mutually absolutely continuous with 
respect to one another; see also Dalang and Nualart ([2004, Corollary 5.3) and Dalang, 
Khoshnevisan, and Nualart (120071 Equations (5.2) and (5.3)). A consequence of this 
mutual absolute continuity is that H b is [Holder] continuous iff H is. 

(2) Theorem 18.21 implies facts that are cannot be described by change-of-measure meth- 
ods. For instance, it has the striking consequence that with probability one, H b and 
H blow up in exactly the same pointsl [This is simply so, because H b — H is locally 
bounded.] For an example, we mention that the operators considered Example 15.51 
when the parameter a there is < 2, lead to discontinuous solutions H that blow up 



a.s.) in every open subset of R + x R ( |Marcus and Rosen, 2006[ Section 5.3). In 



those cases, H b inherits this property as well. 
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Proof. We will need the following fact: 

(8.8) H is continuous in probability. 

In fact, we prove that H is continuous in L 2 (P). 
Owing to (Oil , for all t > and x, y G R, 

,8,) BW,,,-^^^^^* 

Moreover, according to (I3.17p . (1 + Re \l/) _1 G L 1 (R). Therefore, the dominated convergence 
theorem implies that for all T > and t/GR, 

(8.10) sup \\H(t , x) — H(t , y)\\ L 2(p) — > as x — > y. 

te[o,T] 



1 f 00 

(8.11) \imsuv\\H(t,x)-H(s,x)\\l 2{P) <]im— / — 



Similarly, Theorem 15.31 implies that for all £ > 0, 

d£ 

\t-s\) + Re^^y 

which is zero by the dominated convergence theorem. This and ( 18.101) together imply (18. 8p . 
Now we begin the proof in earnest. 
Throughout, we fix 

(8.12) Lip^sup 1 ^' 6 ^ 1 and A := 2Lip fe . 

x ^y \x y\ 

The condition that b is globally Lipschitz tells precisely that Lip 6 and/or A are finite. 

Now we begin with a fixed-point scheme: Set Uo(t,x) := 0, and define, iteratively for all 
integers n > 0, 

i*t poo 

(8.13) u n+l (t ,x) := H(t ,x) + / / b(u n {s ,x - y))p t _ s {y) dyds. 

JO J-oo 

Consider the processes 

POO 

(8.14) D n+ x(x) :— / e~ xt \u n+1 (t,x) - u n (t ,x)\ dt for all n > and x G R. 

Jo 

Also, define r\ to be the A-potential density of X, given by 

POO 

(8.15) r x (z):= p s (z)e~ Xs ds for all z G R. 



According to Lemma 18.11 this is well defined. 
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The processes Di, D 2 , . . . satisfy the following recursion: 

D n+1 (x)<^(D n *r x )(x) 

(8.16) X A 

as can be seen by directly manipulating (|8.13|) . We can iterate this to its natural end, and 
deduce that 

(8.17) D n+1 (x) < 2- n - 1 (D *r x )(x) for all n > 0. 

Since u\{t , x) — u (t , x) = H(t , x) + tb(0), 

r°° w \h(ri\\ 

(8.18) D Q (x)< I e- xt \H{t,x)\dt + 



A 

Thanks to (18.81) . we can always select a measurable modification of (u , t , x) i— > H(t , x)(u>), 
by the separability theory of Doob (119531 Theorem 2.6, p. 61). Therefore, we can apply the 
preceding to a Lebesgue-measurable modification of t i— > H(t , x) to avoid technical problems 
( Doob, 1953 , Theorem 2.7, p. 62). In addition, it follows from this and convexity that 

-oo poo 



/OO POO POO 

\\Do(x)\\mp) e ~ lxldx < / e ~ AHx ' \\H(t,x)\\ L 2 (P) dtdx + 
-oo J — oo J 



A 



We can apply Proposition 13. 1[ with its A replaced by (4A) 1 here, to find that 

\\H(tx)f < eXt/2 r di 

(8.20) \m^)W^< 4?r i_ oo (4A)-i + Re^(0 

:= const • e xt/2 , 

where the constant depends neither on t nor on n. Consequently, 

POO 

(8.21) / \\D (x)\\ L 2 {P) e' lxl dx < oo. 



Because Xr x is a probability density on R, it follows that for all integers n > 0, 

/oo 
\\ D n+i(x)\\ L2{p) e'^dx < const ■ T n . 
-oo 

Therefore, J2^=o ll-Cn(^c) ||z, 2 (P) exp( — [a;|) dx < oo. Furthermore, for all T, k > 0, 

pT pk pk 

lim / / \u n+ i(t ,x) - u n (t,x)\ dxdt < e XT+k lim / D n (x)e~^ dx 
n ^°° Jo J-k n ^°° J-k 



(8.23) 



a.s. 
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Because L x ([0 , T] x [—A; ,k]) is complete, and since (Q , JP, P) is complete, standard arguments 
show that there exists a process Uoo G L] 0C (R + x R) such that 

pT rk 

(8.24) lim / / \u n (t ,x) — u^t ,x)\ dxdt — 0, 

n ^°°Jo J-k 

almost surely for all T, k > 0. Since b is globally Lipschitz, it follows easily from this that 
outside a single set of P-measure zero, 

i*t poo 

(8.25) Uoo(t,x)=H(t,x)+ / b(u 00 (s,x-y))p t -s(y)dyds, 

JO J-oo 

simultaneously for almost all (t,x) G R + x R. An application of Fubini's theorem implies 
then that the assertion (18.25P holds a.s. for almost all (t ,x) G R+ x R. [Observe the order 
of the quantifiers.] 

Consider the finite Borel measure 

(8.26) T(dt dx) := e ~ AHa: l dt dx, 

defined on R + x R. Our proof, so far, contains the fact that G L X (T) almost surely. 
Moreover, thanks to (I8.25p . if p G [1 , oo) then 

(8.27) || Woo j x ) IIlp(p) , x) ||xj>(p) "I" SU P \b{z) \t for all (t,x) G R+ x R. 

z 

One of the basic properties of centered Gaussian random variables is that their pth moment 
is proportional to their second moment. Consequently, 



(8.28) \\ Uoo (t 

i x ) IIlp(P) — cons t ■ ( ||-^(^ ) x ) \\l,2(p) ~l~ 

which we know to be locally bounded. This and the Tonelli theorem together prove that 

Moo e ^L( R + x R ) a - s - 

We recall a standard fact from classical analysis: If f G L 1 ^), then f is continuous in 
the measure T. This means that for all 5 > 0, 

(8.29) limT{(t,x): \f(t + r) , x + e) - f(t , x) \ > 8} = 0, 

and follows immediately from standard approximation arguments; see the original classic 
book of Zygmund (119351 §2.201, p. 17), for instance. Consequently, Lemma I8TT1 and the 
already-proved fact that G L l (T) a.s. together imply that is continuous in measure 
a.s. Because b is bounded and Lipschitz, the integrability /continuity properties of Pt(x), as 
explained in Lemma [87T] together imply that (t , x) \— > f Q 6(^00(5 , x — y))pt- s {y) dyds = 



L b^oci^s , z)) Pt- S (x — z) dz ds is a.s. continuous, and bounded on [0 , T] x R, for every 
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nonrandom and fixed T > 0. Now let us define 

/>t /*co 

(8.30) H b (t ,x) := H(t ,x) + / 6(«oo(s ,x - y))p t - s (y) dyds. 

JO J -co 

Thanks to flQj) , 

(8.31) P{i2" 6 (t,ar) = a:)} = 1 for all (t , ar) G R+ x R. 

Thus, is a modification of Mqo. In addition, the Tonelli theorem applies to tell us that Hf> 
inherits the almost-sure local integrability property of Uoq. That is, Hb G L p loc (Tl + x R) a.s. 
Moreover, outside a single null set we have 

ft /"CO 

(8.32) J(t,x) = / / ^(^(s ,x — y))p t - s (y) dyds for all (t , x) G R + x R. 

J0 J -00 

This proves the theorem. □ 
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